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1 Introduction 



Renewed interest in exploiting the non-perturbative properties of gauge theories with extended 
supersymmetry has been triggered by recent work of Seiberg and Witten [1] (for earlier refer- 
ences see these papers, and for generalizations, see [2]). In order to understand the implications 
of these beautiful, exact results for lattice gauge theory, the lattice regularization of theories 
with extended supersymmetry has to be investigated. This can be done, as proposed some time 
ago by Curci and Veneziano [3], by letting the lattice spoil supersymmetry in the cut-off theory 
and requiring its restoration in the continuum limit. In the case of N = 2 and N = 4 extended 
supersymmetry the definition of the lattice theory is straightforward, because the difficult prob- 
lem of reconciling the lattice with chiral symmetry does not arise due to the vector-like nature 
of these theories. (For a recent review of the chiral problem, see [4].) 

The simplest prototype model with extended supersymmetry is the N = 2 supersym- 
metrized version of Yang- Mills theory with the SU(2) gauge group (SYM2). On the lattice this 
is embedded in a gauge model containing, besides the triplet gluinos, also scalar triplet fields. 
These latter have Yukawa couplings as well as quartic scalar couplings [5] ; therefore this model 
belongs to the large class of Higgs- Yukawa models enjoying interest in the non-perturbative 
lattice literature (for a review see, for instance, [6]). The most general renormalizable quantum 
field theory containing SYM2 has seven couplings and three mass parameters. This makes the 
problem of parameter tuning for the continuum limit quite hard, even in this simplest example. 
Nevertheless, in the weak coupling region a very valuable guide is provided by the perturbative 
renormalization group equations (RGEs). Another useful tool in weak coupling theories with 
scalar fields is the perturbative effective potential, which gives information on the phase struc- 
ture. The usefulness of lattice perturbation theory at weak couplings in lattice Higgs models is 
well known and widely explored (for a recent example, see [7]). In the present paper one-loop 
lattice perturbation theory will be used for the study of the SU(2) adjoint Higgs- Yukawa model 
in the vicinity of the supersymmetric fixed point. 

Among the seven couplings of the renormalizable extension of SYM2, there are four quartic 
couplings. This extended set of quartic couplings plays a dual role: First, it is a complete 
set of couplings which is generated by quantum corrections from any smaller set of couplings 
once supersymmetry is broken. Second, the supersymmetric combination of couplings has flat 
directions which have to be cut-off in order that the path integral be well-defined. 

An important condition for the existence of a supersymmetric continuum limit is the exis- 
tence of a fixed point of coupling ratios corresponding to the supersymmetric relations. It is 
also necessary that this supersymmetric fixed point can be reached from at least one attractive 
direction in the ultraviolet. It has been shown in ref. [5] that both these conditions are fulfilled 
and there is exactly one attractive direction in the region with stable path integral. A closer 
investigation of the one-loop RGEs shows that the supersymmetric fixed point is accompanied 
by another fixed point of coupling ratios, which is actually an infrared fixed point. This will 
be shown in the next section of this paper. 

The Feynman rules of lattice perturbation theory are defined in section 3. Section 4 is 
devoted to the derivation and numerical study of the one-loop lattice effective potential. In 
section 5 questions related to the restoration of global symmetries are investigated. The last 
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section contains a discussion and the conclusions. 

Throughout this paper the same notations will be used as in ref. [5] and some of the equations 
of this reference will only be referred to, without repeating them here. 



2 Infrared fixed point of coupling ratios 

The full set of one-loop RGEs for coupling ratios has been given in eq. (23) of ref. [5]. The 
coupling ratios are defined as 

r, Ga,B ^A,B \aB] \aB) 

Ra,b = , r A>B = — =- , r [AB ] = — =- , r [AB ) = — =— . (1) 

9 9 9 9 

Here g is the bare gauge coupling, G a ,b are Yukawa couplings and X a ,b,[ab],(ab) quartic scalar 
couplings (for the lattice action see next section). The one-loop RGEs preserve the relations 

R A = R B = 1 , r A = r B = r . (2) 

Furthermore, if these relations are fulfilled, we have 

j t [±r A r B - (r [AB] - r {AB) f] = -| [Ar 2 - (r [AB] - r {AB) f] 

= 704r 3 o - r [ll2(r [AB] - r [AB) f + 8 - 327-^J + A(r [AB] - r [AB) ) 

-32(r [AB] - r (AB) ) 3 - lQr 2 (AB) (r [AB] - r [AB) ) . (3) 
The right-hand side here vanishes for 

f[AB] = r (AB) + 2r . (4) 

This implies that relation (4) is also preserved by the one-loop RGEs. It also gives a convergent 
path integral for positive r and non-negative r[ AB y 

If the conditions in (2) and (4) are imposed, the one-loop RGEs are reduced to the simple 
form 

dg 2 _ 8g 4 
dt 16tt 2 ' 

dr g 2 



f 

dt 16?r 
The solution of the first equation is 



^ = 4* ( Ur UB) + Mr r (AB) - 3) . (5) 



a 2(t) = 9 ^ to) (6) 

9{) l + (t-t )g*(t )/(2n*)- {) 
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The last two equations can be written in the variable r = log(<7o/<7 2 ) as 

77 = I ( 112r o + Hab) + Wr r (AB) - l) , 



dr 8 



dr 8 

The right-hand side of the equations in (7) vanish in the points 

S ■ r = , r {AB) = ^ ; (8) 

and 

Vl05 7105 
Q : r = = 0.048795... , r {AB) = = 0.341565... . (9) 

There are two other solutions obtained by reflections about the origin in the (r , r^^-plane, 
but they do not correspond to convergent path integrals. The fixed point in (8) is the super- 
symmetric one in the special case of eqs. (2), (4). Note that the above construction implies 
that every fixed point of (7) supplemented by (2) and (4) is at the same time also a fixed point 
of the full set of RGEs in ref. [5]. 

The stability properties of the fixed points (8) and (9) can be deduced from the derivative 
matrix of the right-hand sides in (7): 



D 



/ 28r + 2r^ AB ) r^ AB ) + 2r 

(10) 



The eigenvalues cti i2 and eigenvectors ei ]2 of this matrix in the supersymmetric point (8) are 
the following: 

a 1 = -l/2 , e 1 = (l,-3), 

a 2 = 3 , e 2 = (l,4) . (11) 

Since, according to eq. (6), for t — > oo also r — > oo, in the continuum limit the negative 
eigenvalues are attractive and the positive ones repulsive. The direction ei corresponds to the 
attractive direction in the supersymmetric fixed point found in ref. [5]. In the other fixed point, 
given by (9), we get 

«i = (5^105 + 3 Vl61)/28 = 3.1893... , e x = (1, -1/2 + ^345/6 = 2.5957...) , 

a 2 = (5^105 - 3 v / 161)/28 = 0.4703... , e 2 = (1, -1/2 - ^345/6 = -3.5957...) . (12) 

These are both repulsive in the ultraviolet and hence attractive in the infrared. The existence 
of this infrared fixed point in the RGEs of coupling ratios has interesting consequences for the 
supersymmetric continuum limit. These will be discussed in section 6 in connection with the 
phase structure at weak couplings. 

The two fixed points S in (8) and Q in (9) determine the flow pattern of renormalization 
group trajectories in the physically interesting range r > 0, T(ab) ^ — 2r . For r ,r^ AB ^ > this 
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Figure 1: The renormalization group trajectories in the (r , r^^-plane around 
the two fixed points S and Q. The latter is attractive in the infrared limit. This 
determines the flow directions. 
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is illustrated by numerical solutions of eq. (7) in fig. 1. The figure nicely shows the separating 
line of the flow which connects the infrared fixed point Q with the supersymmetric fixed point 
S. This line leads to the continuum limit with supersymmetric coupling ratios. It approaches 
S from the ultraviolet attractive direction ei in (11). The flow line approaching S from the 
infrared attractive direction e 2 is also shown. It separates the flow lines going in the infrared 
to Q from those going to the negative r region. 



3 Lattice perturbation theory 

The physical fields in SYM2 are: A r Xfl , fc, W x , K- The gauge field A r Xfl , p G {1,2,3,4}, 
r G {1,2,3} is represented on the lattice by the SU(2) matrix on links U XI1 = exp{igT r A x ^), 
with the SU(2) generators T r = r r /2. The lattice spacing (usually denoted by a) is set to 
1 throughout this paper; in other words every dimensional quantity is measured in lattice 
units. The fermion fields can be represented by the four-component Dirac field in the adjoint 
(triplet) representation ip x , ip x , r, s G {1,2,3}. The real scalar and pseudoscalar triplet fields 
are denoted by A r x , r G {1,2,3} and B s x , s G {1,2,3}, respectively. Given this set of fields, 
we shall consider the most general renormalizable quantum field theory, which respects gauge 
invariance and certain global symmetries, as e. g. parity conservation. This is the SU(2) 
Higgs- Yukawa model studied in ref. [5]. 

The lattice action has been given in eqs. (l)-(4) of ref. [5]. It can be written as 

S = S g + S f + S s . (13) 
The standard Wilson action for the gauge field S g is a sum over the plaquettes 

s 9 = PY.( l -\^ u v) , (14) 

with the bare gauge coupling given by (3 = 4/g 2 . In a triplet notation and with the usual 
perturbative field normalization, the fermionic part Sf is 

Sf = E { ("V + ±r)Tj>l + ierstW x (G A A' x + iG Bl ,B s x )^ x 

X 

~\ E [W x+ »Vrs, x »(r + 7„ty: + W x V-^(r - 7 „)^ + J } . (15) 

Here is the bare fermion (Dirac) mass, < r < 1 is the irrelevant Wilson parameter 
removing the fermion doublers in the continuum limit, and the 0(3) matrix for the gauge-field 
link is defined as 

V rs , x , = 2Tv(U x J r U x M = v; w = V~Z . (16) 
The scalar part of the lattice action S s is 

s s = E { k m A + 8)a:a: + \{m\ + 8)b:b: - £ [ai +p y„^ + Bi + ~v rs , x »B x ] 

x z z n=l 
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+xa(a:ad 2 + \ b (b:b:) 2 + x [ab] a:a:b x b x - \ ab) (aibd 2 } . (17) 

Here m A and m B denote the bare mass parameters for the two scalar fields, respectively. 

In the symmetric phase with vanishing vacuum expectation values the propagators for the 
scalar, fermion and gauge fields are, respectively, 

XUX2 (m^ + fP^ + P V ^2 ^ 7 j ' 

, 1 r i i 1 

(18) 



AS _ ^ Jk-(xi-x 2 )+i(k u , -k uo )/2 



<W 2 ~ I 1 ~ a ) 



where, as usual, f2 is the number of lattice points, a is the gauge-fixing parameter and 

^ = 2sin(^), ^ = sin^) . (19) 



The vertices can be read from the action (13)— (15) after expanding the gauge links, with 
\Axn\ 2 = A r Xf ,A r x ^ according to 

a: 



U XIX = cos (Jt + 2zT r |-^| sin (l^l) = a° x/x + 2iT r a r x/x ; 



XfJ, 

2 



= S r . + ge TSt A^ + 9 - (Al^ - £, S <X J + . . . . (20) 

The Feynman rules in the phases with non- vanishing vacuum expectation values 

v A S Sr = (Al) + , v B S Sr = (Bl) + (21) 

are, of course, somewhat more complicated. As has been discussed in ref. [5], according to 
the classical potential there are three different phases with the symmetry-breaking pattern 
SU(2)— s-U(l): phase A with v A 7^ 0, v B = 0, phase B with v B ^ 0, v A = and phase AB with 
v A 7^ 0, v B ^ and parallel vacuum expectation values as in (21). (The phase with orthogonal 
vacuum expectation values and complete breaking of the SU(2) gauge symmetry is not relevant 
to supersymmetry.) It can be expected that for weak couplings the phase structure remains 
qualitatively the same. In fact, this is the main question which will be investigated in the 
present paper within the framework of lattice perturbation theory. 

In order to illustrate the Feynman rules in these phases with Higgs mechanism, let us briefly 
consider the most interesting case of phase AB. The other cases are simpler and can be treated 
similarly. Here only a few remarks will be included about them, which are relevant in some 
important, degenerate situations. Let us split the scalar field according to 

Ai = A: + v A 8 3r , b: = b: + v b 8 3i , (22) 
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and consider the gauge boson propagator. The gauge-fixing function is conveniently introduced 
in such a way that the scalar-gauge mixing is cancelled [8]. This is achieved by taking, for the 
gauge-fixing part of the action, 

Sgf = frxfrx (23) 



2a 



with the gauge-fixing function 



f rx = ^2(A r xli - A r x _~ ^) + age rU {v A A x + v B B x ) 
With this choice the gauge boson propagator becomes 



(24) 



A 9 =-Ve 



'.k-(xi —X2)+i{k lll — fe M2 )/2_ 



h h 

W 2 - I 1 - a ) — : — -77 

am g,r + ^ 



Here the gauge boson masses m giT are given by 

2 2 2/ 2 1 2 \ 

V = m g,2 = 9 \ V A + V B) 



<3 = 



(25) 



(26) 



The mass square matrix for the scalar fields (A, B) can be obtained by introducing the 
shifted fields (A, B) defined in (22) into the lattice action and collecting the quadratic terms 
in the minimum of the potential. Adding also the contribution of the gauge-fixing term pro- 
portional to a, one obtains for the mass square matrix in the subspace with isospin indices 1,2 
containing the Goldstone bosons: 



M 2 



( 2\ab)V b + ag 2 v\ -2\ AB )V A v B + ag 2 v A v B 

(12) = 

V ~2\ab)Vav b + ag 2 v A v B 2\ab)V a + ag 2 v 2 B 



(27) 



(28) 



In the orthogonal subspace with isospin index 3 we get 

( &\ A v 2 A ±{\ab] - \ AB ))v A v B 

V 4:(\ab] - \ab))vavb 8A b v| 
In the Landau gauge (a = 0), these give for the momentum space scalar propagator 

G AB = (M ( 2 123) + k 2 )- 1 (29) 

the non- vanishing matrix elements 

fiAB _ nAB _ A; 2 + 2\ AB) v 2 A 

^11 - ^22 



k 2 [k 2 + 2\ {AB) (v 2 A + v 2 B )] 



^AB _ ^AB _ k 2 + 2\ AB) V 2 B 



k 2 [k 2 + 2\ AB) (v 2 A + v 2 B )] 



f,AB _ f,AB _ f,AB _ f,AB _ 1\ AB )V A V B 
Lt 14 — Lt 25 — Lt 41 — Lt 52 



k 2 [k 2 + 2\ {AB) {v 2 A + v 2 B )} ' 
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fiAB _ k + 8X BV 2 B 

^33 



G 6 6 



k 4 + 8k 2 (X A v 2 A + X B v 2 B ) + Wv 2 A v 2 B [AX A X B - (X [AB] - X (AB) y 
ab _ k 2 + SX A v A 



k 4 + 8k 2 {X A v 2 A + X B v 2 B ) + Uv A v%[AX A X B - {X [AB] - X (AB) y 



qAB = qAB = -*(\ab] - \ab))vavb , 3Q , 

63 k 4 + 8k 2 (X A v 2 A + X B v 2 B ) + 16v 2 A v 2 B [AX A X B -(X [AB] -X (AB) ) 2 ] ' 

The indices are denned here in such a way that the first three belong to A and the last three 
to B. The notation k 2n means, of course, [k 2 ) n . 

An important, special, degenerate case of the phase structure occurs when the All-phase 
collapses to a line in the (m 2 A , ra^)-space and, therefore, the A- and 5-phases touch each other. 
At tree level this happens when 

2 I \ 2 



^X A X B = (\[ab] - \ab)) 5 m A \JX B = m B \JX A . (31) 

In this case, on the line unifying all three phases, there is a minimum valley where the vacuum 
expectation values of the A- and .B-fields are on an ellipse 

v A <J X A + v B <J X B = - =- . (62) 



WX A 4vA; 



In fact, the perturbation theory in this degenerate situation can be considered as a special 
case of the perturbation theories in any of the three phases. For defmiteness, one can transform 
the vacuum expectation value to the yl-direction by an axial U(l)yt transformation 
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A' x = cos(2(j))A x - sm(2(j))B x , B' x = sm(2<j>)A x + cos(2<j>)B x . (33) 
Then, for instance, (30) becomes 

fiAB _ fiAB _ fiAB _ 1 fiAB _ \ fiAB _ fiAB _ \ (o,\ 

k 2 k 2 + 8X A v A k 2 + 2\( AB )V A 

As one can see, besides the Goldstone bosons A 1 , A 2 there is also a third zero-mass boson B 3 . 
This is one of the scalar N = 2 superpartners of the massless photon field . Since it is 
associated with the global axial U(l)yt transformation, we can call it photo- axion. The other 
scalar superpartner of the photon field (A 3 ) is an excitation orthogonal to the minimum valley 
in the plane of the two vacuum expectation values. Since it has to do with the length of the 
vacuum expectation values, we can call it photo- dilaton. In the same spirit, the scalar members 
of an N = 2 vector (gauge-) supermultiplet can be called gauges cellars, gluo-scalars or photo- 
scalars. In perturbation theory the photo-dilaton is not exactly massless, but as can be seen 
from (34), its mass square is proportional to A^. Near the supersymmetric continuum limit 
g 2 — > 0, according to eqs. (7)— (11), goes to zero as A^ oc g. Of course, these are just masses 
in the propagators of perturbation theory and the real question is how the physical masses 
of the photo-axion and photo-dilaton behave in the continuum limit near the supersymmetric 
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fixed point. For answering this question, one has to first see whether the appropriate phase 
structure is reproduced at all near the continuum limit (see next section). 

Returning to the Feynman rules, in the AB -phase the inverse propagator of the fermion 
field in momentum space is 



Gt.l = 8 rir2 ( m^p + -k 2 + ij ■ k ) + ie ri3r2 (G A v A + ijsGBVs) 



(35) 



After introducing the notations 



Gf 



m^ k + G a va - ij5G B VB - i"1 ■ k 
(m^u + G A v A ) 2 + G 2 B v 2 B + k 2 

_ m^k - Z7 -k 



Gl 



m^ k - G A v A + ij 5 G B VB - i"1 ■ k 
(ra^ - G A v A ) 2 + G\v\ + k 2 



Gl 



m^,k = + -k 2 , 



(36) 



the fermion propagator in momentum space can be written as 

/ \{Ot + Ot) UGt-Gt) \ 



G^ 



HGt + Gi) o 



v 











(37) 



Gt J 



Besides the propagators, one also has to determine the vertices. They can be obtained also 
in the Higgs phases by simply expanding the lattice action in powers of the fields A r in (20), 
A;, 5; in (22), and ^J x . 



4 Effective potential and phase structure 

The existence of a supersymmetric continuum limit in a renormalizable lattice theory, with the 
same set of fields and a set of couplings and mass parameters broader than the supersymmetric 
target theory, imposes several conditions on the behaviour of the bare theory for small lattice 
spacings. In the case of an asymptotically free continuum limit, as required for SYM2, the 
continuum limit is defined by approaching the free Gaussian fixed point at zero couplings. The 
first condition to be satisfied is the existence of a fixed point in the RGEs for coupling ratios 
corresponding to the supersymmetry relation between the gauge coupling (g) and the other 
relevant couplings. This fixed point should occur for zero gauge coupling and should have at 
least one attractive direction in the g — > limit, in order that it can be reached from non-zero 
g. As has been shown in ref. [5], this condition is satisfied in SYM2 (see also section 2 of this 
paper). 

The second condition for the existence of the supersymmetric continuum limit is that the 
appropriate phase structure, which arises in the classical theory from the classical potential, 
is also reproduced near the Gaussian fixed point in the quantum theory. Besides the trivial 
requirement of second-order nature of the phase transitions, which makes the continuum limit 
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possible in general, there are also more subtle conditions to satisfy. These follow from the ne- 
cessity of reproducing the flat directions in the effective potential corresponding to the quantum 
moduli space (QMS) studied in ref. [1]. The classical moduli space (CMS) of SYM2 arises as a 
particular limit of the minima of the classical potential in the SU(2) Higgs- Yukawa model. The 
compact flat direction belonging to the angle of the vacuum expectation values of the scalar 
A- and 5-fields corresponds to the minimum valley, which appears if the phases A, B and 
AB become degenerate on a line in the plane of bare scalar mass squares (m^,m^) (or scalar 
hopping parameters {^a-,^b))- The non-compact flat direction belonging to the overall length 
of the vacuum expectation values is obtained in the limit where the quartic scalar couplings 
stabilizing the path integral (A^, A#) tend to zero. As already discussed in the previous section, 
the existence of these two flat directions in the QMS requires that the physical mass of the 
photo-axion and photo-dilaton be zero in the continuum limit. In summary: the requirement on 
the phase structure is the existence of second-order phase transitions having the same topology 
as in the classical theory, in particular, showing the degeneracy of the three phases A, B and 
AB in a line of the (m A , ra^)-plane. 

The phase structure of lattice quantum field theories can be effectively exploited by nu- 
merical simulations. Another useful tool, which is reliable in models with scalar fields with 
weak couplings, is the study of the effective potential in lattice perturbation theory. Previous 
experience shows that the nature and position of phase transitions in bare parameter space is 
quite well reproduced by the one-loop approximation. For instance, in the SU(2) fundamental 
Higgs model in the physically relevant region (3 ~ 8 - 10 and A ~ 10~ 4 - 10~ 3 both the invariant 
effective potential (in the unitary gauge) and the Landau gauge (a = 0) effective potential work 
well [7]. (In fact, in the published version of ref. [7] only the invariant effective potential is 
considered, but very similar results can be obtained in the Landau gauge too.) In the present 
context the study of lattice perturbation theory is even more important than in simple Higgs 
models. This is because the numerical simulation is much more expensive, due to the presence 
of fermions and the large number of parameters. In the present section I shall consider the 
one-loop lattice effective potential in the Landau gauge. 

The one-loop effective action for the bosonic fields is given in general by 

r^op^,^,^] = S b [Al^Al,B:} + ±TT\og{D b [A:^Al,B:}A b } 

-TT\og{M FP [A:^Al,B:}A FP } - TT\og{M*[A:^A:,B:}A*} . (38) 

Here Sb = S g + S s is the bosonic part of the lattice action in eqs. (13)— (17), Db its second 
derivative matrix, A 6 the bosonic propagator matrix, M FP the Fadeev-Popov matrix for the 
ghost fields with the ghost propagator matrix A FP , A^ the fermion propagator matrix and 
the fermion matrix in the action defined by writing the fermionic part as 

s>= E 7yK.r.K ■ (39) 

rx,sy 

The one-loop effective potential is defined for ^-independent scalar fields A r ,B r by 

Kff" loop = H r i-wK, = o> A : = Ar , b : = Br ] ■ (40) 
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Up to now these are the unrenormalized expressions. In the vicinity of second-order phase 
transitions, where the masses and vacuum expectation values tend to zero in lattice units, 
the bare perturbation theory in not convergent. This is due to the large logarithms log(ara), 
where am is a typical mass or vacuum expectation value in lattice units. Here we would like 
to determine the phase structure in bare parameter space; therefore, we shall not completely 
switch to renormalized perturbation theory, but only replace the masses in the propagators by 
the renormalized ones. This can be achieved in the loop expansion by adding to the action 
appropriate external sources proportional to the mass terms, which are quadratic in the bosonic 
fields. It is a general experience in lattice quantum field theory that the location and nature of 
phase transitions in bare parameter space are then well reproduced. In fact, an even simpler 
approximation is working well, when the masses in lattice units are set to zero everywhere 
in the propagators. This is the approximation I shall exploit in what follows. The effective 
potential will be calculated in the symmetric phase without scalar field expectation values and 
will be studied numerically in the vicinity of the phase transition hypersurfaces. 

The technical advantage of considering the Landau gauge, with gauge parameter a = 0, 
is that there is no mixed scalar-gauge loop contribution and the Fadeev-Popov ghost loop 
contribution vanishes, as well. Hence the one-loop effective potential in the Landau gauge can 
be written as the sum of scalar, fermion and gauge loops: 

x^l-loop T/scalar , T/fermion , T/gauge / . -. \ 

K eff — K eff ' V eS T ^eff ■ 



Let us first consider the scalar contribution, which is in fact the most complicated one out 
of the three. According to eqs. (38)-(40) it can be written as 

=2?T? l0gdet i 1+ (^ V BB )[ (mi + ^-'Jl 

^EloglV (42) 
It follows from the lattice action in eq. (17) that the 3 ® 3 coupling matrices V. are given by 
V AA , rir2 = 8 rir2 (a\ a {AA) + 2\ AB] {BB)) + &\ A A r 'A r > - 2\ {AB) B ri B r 

VBB, ri r 2 = S T1T2 (A\ B {BB) + 2\ [AB] {AA)) + 8\ B B ri B r2 - 2\ AB) A ri A r 

V AB>riT2 = 4:X[ AB ]A ri B T2 — 2\ AB )8 riT2 (AB) — 2\ AB )B ri A 7 " 2 , 

V BA<rir2 = AX [AB] B ri A r2 - 2\ AB) 8 rir2 {AB) - 2\ {AB) A ri B r2 . (43) 

Here for the scalar products of isospin vectors the notation (AB) = A S B S etc. has been intro- 
duced. The evaluation of the determinant T> s in eq. (42) is possible by an algebraic manipulation 
program, but the resulting expression is rather long. In the important special case \ B = \ A , 
\ab] = A( J 4B)+2A J 4, which corresponds to the relations (2), (4) for the coupling ratios considered 
in section 2, the result is given in the appendix. Further specializing to 

g 2 

As = \ A = , \ A b] = \ab) = — , (44) 



?r 2 



\r 2 

5 
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we get 

V s = 1 + 2g 2 D 1 + g\D\ -W 2 )- Ag 6 D 1 D 2 - g\D\ + D 2 )D 2 + 2g w D 1 D 2 2 + 3g 12 D 3 2 . (45) 
Here the abbreviations 

Dl = Dl(A , B , k) =M + (BB) 



m 2 B + k 2 m 2 A + k 2 

, N (AA)(BB) - (AB) 2 , . 

D 2 = D 2 (A, B, k) = 2 *;> 2 \^ (46) 
(m A + k 2 )(m B + P) 

have been used. 

The fermion-loop contribution in eq. (41) can be similarly calculated to give 

K f |rmi on = _2 £ logI)/ ^ (47) 

where 

V f = l + D} [2G 2 A (AA)(k 2 - mj fc ) + 2G 2 ,( J B£)(P + mj fc )] 

+0} [G}(AA) 2 (* 2 + m^) 2 + G%(BB) 2 (k 2 + m\ k ) 2 

+2G 2 A G 2 B (AA)(BB)(k* - mj fc ) + 4G 2 4 G 2 ,(^) 2 (P + m 2 fc )m 2 J (48) 
and with (36) we defined 

^ = (m 2 fe + P)" 1 . (49) 
Finally, the contribution of the gauge-boson loop is given by 



where 

Here the notations are: 



V g = 1 + 2g 2 D' 1 + 5 4 (Z)f + £> 2 ) + g*D' x D' 2 . (51) 

(AA) + (BB) 



Jfe 



2 



(52) 

The one-loop expression in eq. (41) for the effective potential in the Landau gauge can 
be evaluated numerically for any choice of the couplings and invariant scalar field variables 
(AA), (BB) and (AB) 2 . As discussed before, for the study of the phase structure at weak 
bare couplings a good approximation can be expected, if one replaces the bare masses in the 
propagators by the renormalized ones. As a first approximation, one can further simplify the 
analysis by setting the propagator masses to zero: tua = rns = = 0. (Later on this can 
be refined, for instance by combining numerical information on the renormalized masses with 
these perturbative formulas.) The numerical work can still become voluminous even in this 
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Table 1: The parameters of the points where the numerical study 
of the one-loop effective potential was performed. 



Label 


g 2 


^0 


r {AB ) 


\a 


\ab) 


a 


0.5 


0.02 


0.44 


0.01 


0.88 


b 


0.5 


a/105/210 


a/105/30 


0.02440.. 


0.1708.. 


c 


0.4 


a/105/210 


a/105/30 


0.01952.. 


0.1366.. 


d 


0.2 


0.02 


0.44 


0.004 


0.088 


e 


0.2 


0.05 


0.35 


0.01 


0.07 



simplified case, if large lattices are considered, which is important in order to diminish finite 
volume effects in the position and nature of the phase transitions. 

The numerical calculations were performed on 32 4 and 64 4 lattices, after the initial search 
on smaller lattices. The coupling ratios were always chosen to satisfy the relations in (2) and 
(4). The gauge coupling was in the range 8 < (3 < 20 (that is 0.5 > g 2 > 0.2), in order to stay 
within the weak coupling region. The parameters of five points, where most of the calculations 
were performed, are collected in table 1. The values of the scalar-mass parameters m\ B in the 
tree level potential 

Kff ee = ^(AA) + ^(BB) + \ A {AAf + \ B {BBf + \ [AB] (AA)(B B) - \ AB) {ABf (53) 

were changed as long as the phase structure was sufficiently clarified. Of course, the bare mass 
squares rn\,rn B can also be given in terms of the scalar hopping parameters k, a ,k, b , which are 
usually preferred in numerical simulations. The relations are: 

ka,b = (8 + mis)- 1 . (54) 

In the case where is set to zero in the fermion one-loop contribution, the effective potential 
does not depend on the fermion hopping parameter n^. Of course, the result is a good approxi- 
mation only for such values of where the renormalized fermion mass in lattice units is small. 
(For a perturbative estimate of the critical fermion hopping parameter, where the renormalized 
fermion mass vanishes, see the next section.) 

For the points specified in table 1 the values of V^ ee + loop were tabulated in an appro- 
priate range of field variables, in order to find the minimum. In the minimal point the values 
of the first and second derivatives were also calculated from the analytic expressions, in order 
to correct for the minimum position and get an estimate of the scalar masses. In fact, choosing 
a fine enough mesh of points for the tabulation, the correction on the minimum position could 
be kept very small. The obtained phase structure in the (k, a , Ke)-plane is shown in figs. 2-6. In 
these figures S denotes the symmetric phase with v A = v B = and A, B, AB the three Higgs 
phases, respectively. 

The dangerous possibility of the appearance of first-order phase transitions is realized at 
the parameter values in fig. 2. In this case the radiative corrections are destroying the clas- 
sical picture completely. Large parts of the renormalized parameter space are cut out by the 
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phase structure at g 2 = 0.5, r= (0.02,0.44) 



0.132 



0.13 



CQ 



0.128 




0.134 



0.136 



0.138 



A 



Figure 2: The phase structure in the plane of scalar hopping parameters (ka,k.b) 
at g 2 = 0.5, r = 0.02, ?"(ab) = 0.44. The dashed line denotes a first-order phase 
transition, the full lines second-order ones. 
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I I I I I I I I I I U I I I I I I 

0.138 0.139 0.14 

Figure 3: The phase structure in the plane of scalar hopping parameters (k.a, Kb) at 
g 2 = 0.5 and the infrared fixed point Q for the coupling ratios. The full lines denote 
second-order phase transitions. 
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I I I I I I I U I I I I I I I I 

0.135 0.136 0.137 

Figure 4: The phase structure in the plane of scalar hopping parameters (ka, Kb) at 
g 2 = 0.4 and the infrared fixed point Q for the coupling ratios. The full lines denote 
second-order phase transitions. 
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I I I 1 I I I I I I U I I I I I 

0.129 0.13 0.131 

/c A 

Figure 5: The phase structure in the plane of scalar hopping parameters (ka,k.b) 
at g 2 = 0.2, r = 0.02, ?"(ab) = 0.44. The full lines denote second-order phase 
transitions. 
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I I I I I U I I I I I I I I I I 

0.13 0.131 0.132 

Figure 6: The phase structure in the plane of scalar hopping parameters (ka,k.b) 
at g 2 = 0.2, r = 0.05, T{ab) = 0.35. The full lines denote second-order phase 
transitions. 
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metastable regions near the first-order transition. Fortunately, in the other cases no first-order 
phase transition is seen. Within the precision of the numerical calculation, no sign of a disconti- 
nuity is seen. All phase transitions are either of second order or at most very weakly first order. 
(Of course, a weak enough discontinuity can never be excluded by a numerical calculation with 
finite precision.) The only qualitative change in figs. 3-6 compared to the tree-level potential 
is the opening up of the degenerate line containing the degenerate AB-phase. (Remember that 
the coupling relation in (31) is satisfied due to eq. (4), therefore at tree level the AB-phase is 
degenerate at all these points.) The conclusion from these figures is that, in order to avoid a 
first-order phase transition, for a given value of the bare gauge coupling g 2 the value of r can- 
not be smaller than some critical value. Otherwise the fluctuations in the lengths of the scalar 
fields become too strong and a first-order phase transition is induced. In the supersymmetric 
continuum limit, however, r — > must be realized. This requires that the critical value of r 
has to become smaller and smaller for g 2 — > 0. Indeed, comparing figs. 2 and 5 shows that 
the value r = 0.02 is well below the critical value for g 2 = 0.5, but comfortably above it at 
g 2 = 0.2. In order to reach a better closing of the AB-phase towards degeneracy at g 2 = 0.2, 
one has to go up to r = 0.05, which is already slightly above the infrared fixed point Q in 
fig. 1. For reaching the flow line connecting the fixed points Q and S with a thin wedge for the 
AB-phase one has to go to an even smaller g 2 . The required tendency of closing the AB-phase 
towards degeneracy for fixed (?"o, ?"(ab)) an( i decreasing g 2 can be seen at Q by comparing fig. 3 
with fig. 4. This is well displayed by the insets of these figures, showing the behaviour near the 
critical point n c where all four phases meet. 

From the point of view of reaching supersymmetry, the best behaviour of the phase structure 
is shown in fig. 6, where the AB-phase is already constrained to a narrow wedge. From the 
experience obtained with moving around in g 2 and (r , 7"(ab))5 it is quite clear that the AB-wedge 
can be closed even further. This has to be checked, of course, later on by numerical simulations. 
Nevertheless, the proof of an exact degeneracy cannot be reached by purely numerical methods. 

Another information that can be obtained from lattice perturbation theory is an estimate 
of the renormalized masses. What is needed is the second derivative matrix M 2 of the ef- 
fective potential with respect to the field components, together with the Z-factors of wave 
function renormalization. The latter were not yet determined. However, general experience in 
Higgs- Yukawa models tells us that the values of these Z-factors are usually not far from 1 [9] . 
Neglecting their deviations from 1, mass estimates can be obtained from M 2 alone. An inter- 
esting question concerns the masses of the photo-axion and photo-dilaton in the AB-wedge, 
which should vanish in physical units in the N = 2 supersymmetric continuum limit. It turns 
out that in situations similar to the inset of fig. 6 the ratio of the photo-dilaton mass to the 
masses of the two massive scalar states is in the range of about 1/2, whereas this ratio for the 
photo-axion is in the range 1/5 to 1/15. As a example, the obtained mass estimates along the 
line k,a + Kb = 0.2590 are shown in fig. 7. 

Both the smallness of the masses and the absence of discontinuities at the phase transition 
lines in figs. 3-6 indicate second-order or very weakly first-order phase transitions. The question 
of the order of the phase transitions needs, however, further study by numerical simulations and 
by refined perturbative analysis. For instance, one could put in the propagators of the pertur- 
bative effective potential some information on the renormalized masses either from numerical 
data or from perturbative calculations at moderately small masses in lattice units. This would 
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Figure 7: The mass estimates obtained from the one-loop effective potential at 
g 2 = 0.2, r = 0.05, T(ab) = 0.35 and k,a + Kb = 0.2590 as a function of k,a- Full 
squares stand for the equal masses of the two massive scalars, open triangles for the 
mass of the photo-dilaton and open circles for the mass of the photo-axion. The 
vertical lines indicate the positions of the phase transitions. 
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improve the ability to disentangle weak first-order phase transitions and real second-order ones. 
The position and nature of the phase transitions has in any case to be investigated by numerical 
simulations, in order to study non-perturbative effects. 



5 Global symmetries 



Besides the vanishing of the masses of the photo-axion and photo-dilaton, another requirement 
in the continuum limit is the restoration of the global SU(2)k symmetry. In the lattice action 
this is broken by the fermion mass term and by the Wilson-term (proportional to r) to U(1)jf 
of fermion number conservation. The SU(2)k symmetry can be seen in the massless continuum 
Euclidean action if one writes it in terms of left-handed fields only. This is achieved by replacing 
the right-handed components, using their relations to the charge-conjugated fields: 

fa = ct cL , = €lG ■ (55) 

Here C is the usual charge- conjugation matrix for Dirac fields, which satisfies 

Cj.C- 1 = - 7 J , C = -C T = -C- 1 . (56) 

The elements U-ji of SU(2)k can be defined by their action on the left-handed fields as 

It is sometimes advantageous to introduce Majorana fields instead of the Dirac fermion field 
used up to now. Their definition is 

= ±=ty + ct) , * (2) = ^(-^ + Ct) ■ (58) 



These satisfy (j = 1,2) 



For the left-handed fields we have 



* 0) = ¥ j)T C . (59) 



where the U(2) matrix S is given by 

*=M-> !) =s ""- (61) 

This implies that, instead of U-ji, the left-handed Weyl-Majorana doublet ty L ' transforms 
under SU(2)k by SUjzS^. Combining eqs. (55)-(60) one can also see that the right-handed 
Weyl-Majorana doublet ^^'^ transforms at the same time by t^SU-jzS^t^, with the Pauli matrix 
r 2 acting on the doublet index (j = 1,2). 
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On the Majorana basis denned by eqs. (58) and (59) the fermionic part of the lattice action 
S7 = E{ krn+ + 4r) £^>¥<?> + ^>{G A K + iG Bl *B' m )*™ 

* Z 3 = 1 

~ \ E E \*%V rs , x »{r + 7 „)*(?> + ^V-^(r - 7 ,)*?> 1 } . (62) 



»=1 j=l 

This form is useful for the derivation of the Ward-Takahashi (WT) identities belonging to 
the SU(2)k symmetry. Performing in the path integral an a;-dependent infinitesimal SU(2)k 
transformation: 

*L = {l-i a ix75Ti+ia 2x T 2 +ia 3x j 5 T 3 )if! x , = ^ x (l-ia lx j 5 T 1 -ia 2x T 2 + ia 3x j 5 T 3 ) , (63) 

with infinitesimal ctj x , the invariance of the path integral measure implies the WT identities. 

First of all, for the unbroken U(1)jf subgroup the gauge-invariant vector current of fermion 
number can be defined as 

= 4? = \ pW>K„(r + 7m)*^ + 9^V-*(r - 7m)**+a] • (64) 
This satisfies the unbroken WT identity 

o = (EK - = (a^) • ( 65 ) 

The other two gauge-invariant currents J^ =1 ' 3 \ which can be defined as 

4? = ^(-i)°'- 1)/2 [*,+a^75^*, + ^wsi;; 1 ^] , (66) 

correspond to the other two generators of SU(2)k. The corresponding symmetries are explicitly 
broken by the fermion mass and by the Wilson term proportional to r; therefore the currents 
41 =1 ' 3 ^ satisfy the broken WT identities 

(AJ4?) = (-l) y " 1)/2 {m^ x r n ^ x + IW) . (67) 
Here the symmetry-breaking term proportional to r is given by 



X i j) = 7 E [4*0,^75*. - V x+i >T jl5 V x ^ x - y x T jl5 Vjy x+il 



-tf^T^sK^tf. - $,r i 7 B K. p ,,* I . jl ] . (68) 

The WT identity in (65) corresponds to exact fermion number conservation at non-zero 
lattice spacing. The broken WT identities in (67) are very similar to the flavour non-singlet 
axial- vector WT identities in lattice QCD with Wilson quarks [10]. Taking into account operator 
mixings, and appropriately redefining the current normalizations and the critical hypersurface 
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for zero fermion mass, one-loop lattice perturbation theory suggests that in the massless contin- 
uum limit these identities are reproducing the expected unbroken WT identities in the target 
continuum theory. 

Of course, the zero fermion-mass hypersurface defined in this way corresponds only to one 
possible definition at non-zero lattice spacings. There are many other definitions, which all 
tend, according to lattice perturbation theory, to the same hypersurface in the continuum 
limit. One possibility is to calculate the fermion self-energy in the AB -phase at vanishing 
four-momentum p = 0. Applying the Feynman rules of section 3 one obtains the general form 

E*(p = 0) = MfX + i^Mf^ . (69) 

The matrices in isospin indices 7"i,r 2 satisfy, at non-zero lattice spacing: 

= mh - mh = mh + mh = mh = mh = mh = mh , 

= Mil - Mil = Mil + Mi 5 = Mil = Mil = Mil = Mil ■ (70) 

In the continuum limit, when in lattice units va,vb — > 0, the only non-zero matrix elements 
are 

TiW'i _ _ 1 ^ 4(^-^-4^-8^-2^ 

11 " 33 " n \ P(4P + jfe*) • (71) 

Here we have set, for simplicity, the Wilson parameter to r = 1. The simplest definition of the 
critical value of the bare fermion mass m c ^ is to require m c ^ + MH = 0. This gives 



m:, = 4 



4(4<7 2 + GI- G\) + 2<7 2 



4Jfe 2 + k A k 2 



(4# 2 + G%- G A ) ■ 0.085.. + g 2 ■ 0.310.. (72) 



Remember that the critical fermion hopping parameter can be obtained from by = 
1/(8 + 2m£). 

Let us note that from one can also extract the one-loop corrections to the renormalized 
Yukawa couplings G T A n B . They turn out to be, respectively, 

lim Mjl _GAy 16g 2 fe 4 + 4(gj - G\ - 4g 2 )(4fe 2 fe 2 - fe 6 ) - g 2 k 8 
vam-+o iv A Q Y k 4 (Ak 2 + 

lim M ^ -° B Y ^ + ^ ~ ^ ~ 4g2) ^ + 4 ^ P (73) 
vam-+o v B Q Y & 4 (4P + A: 4 ) 

The other important global symmetry besides SU(2)k, which has to be restored in the 
continuum limit is, of course, N = 2 supersymmetry. In order to discuss the supersymmetry 
transformations, it is advantageous to introduce the so-called "symplectic Majorana" fields [11], 
which transform simply as a doublet under SU(2)k. These are related to the Majorana fields 
in (58) by _ 

A = {t 2 P r + P L )tf , A ee *(r 2 P L + P R ) . (74) 
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Here Pr = (1 + 75)/2, Pl = (1 — 7s)/2 denote the chiral projectors and the Pauli matrix r 2 
acts on the SU(2)k doublet indices j = 1,2. Instead of eq. (59), these satisfy the symplectic 
Majorana relations 

A = -X T C l5 r 2 , A = C 15 t 2 X T . (75) 

The N = 2 supersymmetric Euclidean action in the continuum can be written with these fields 
as 

5 susy = J d*x{±F;„(x)F;„(x) + l -{D^{x)){D.A\x)) + ±(D ,B' (x))(D ,B' (x)) 

+ \^ jr (^D^ jr (x) - l -^e rst \ jr (x)[A s (x) + i lh B s {x)]V\x) 

+ — [A r (x)A r (x)B s (x)B s (x) - A r (x)B r (x)A s (x)B s (x)]\ , (76) 
where the notations for the field strength and the covariant derivatives are, as usual, 
F; v = d,Al(x) - d v Al(x) + ge rst A^(x)Ai(x) , D,A r (x) = d,A r (x) + ge rst Al{x)A\x) . (77) 

Global N = 2 supersymmetry is realized non-linearly in a formulation without auxiliary 
fields [11]. Introducing the symplectic Majorana doublet Grassmann variable e, which satisfies 
the same relations as (75) 

e = -e T C^ 5 T 2 , e = C^ 5 t 2 I t , (78) 

the supersymmetry transformations are defined as 

6A r = e\ r , 6B r = ie 75 A r , SA^ = ie 7jti A r , 

SX r = - l»{D^A r + i l5 D^B r )e + ge^A'B^e . (79) 

Here = (z/2)^, 7^] is used. By partial integration and neglecting surface terms, it can be 
shown that (79) is a symmetry of the above action: 

£Ssusy = . (80) 

There is a gauge-invariant conserved current corresponding to this symmetry. It is the SU(2)k 
doublet spinor- vector current 

J^(x) = j-I^C^)^ + I^^C^) + i^^C^)]^ - | e _ t ^( a; )^( a; ) T5 }^A-( :c ) . (81) 

As a consequence of the field equations of motion, the divergence of J^ usy vanishes. 

For the study of the supersymmetric WT identities, one has to define the transformations 
(79) on the lattice. This requires some definition of the field strength tensor i^" lce , for instance 
the usual "clover" definition averaging over the four open plaquettes in the /i, i>-plane, which 
are touching the site x. Since the gauge link variables are given by U XI1 = exp(igA r xfJ T r ), the 
transformation of A r can be represented, for instance, by 

SU xli = -ge^YTrU^ . (82) 
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By eq. (16) this also implies the transformation of the gauge links in the adjoint representation 
V rs>xl j,. The covariant derivatives on the lattice are replaced by covariant differences defined, for 
instance, as 

A(U),Al = I [V-^AI + - - K s ,,-a,^:_J • (83) 

All these choices contain, of course, some arbitrariness at the level of O(a) corrections, which 
go to zero by the lattice spacing a. 

Having defined the supersymmetry transformations on the lattice, the broken WT identities 
corresponding to N = 2 supersymmetry can be derived. Due to the involved operator mixings, 
the analysis of the continuum limit is in this case more subtle than in the N = 1 case discussed 
in ref. [3] and will not be considered in the present paper. 

6 Discussion and conclusion 

The flow pattern of the renormalization group trajectories at weak couplings in the SU(2) 
adjoint Higgs- Yukawa model makes it possible to tune the coupling ratios to the values required 
by N = 2 supersymmetry. The fixed-point structure is particularly transparent in the simplified 
version of the renormalization group equations studied in section 2, where several coupling ratios 
are set to satisfy relations that are left invariant by the complete set of one-loop equations in 
ref. [5]. It is interesting that, besides the supersymmetric fixed point of coupling ratios S, there 
is also another fixed point Q, which is attractive in the infrared limit (see fig. 1). This fixed 
point is particularly important for the behaviour of the ratios of renormalized couplings, if the 
bare couplings are kept fixed, because the flow directions are then reversed. In order to avoid 
the attraction of this non-supersymmetric fixed point, one has to reach the separating flow line 
that leads, in the infrared limit, to the fixed point S from the direction e 2 in eq. (11). However, 
it should be kept in mind that in the infrared limit the couplings become strong, and therefore 
the unknown higher-order corrections become important. 

The N = 2 supersymmetric Yang- Mills theory has a weak coupling limit, when the ex- 
pectation values of the gauge-scalar fields get large compared with the A-parameter of the 
asymptotically free gauge coupling. This is a great advantage, because in this region the theory 
is perturbative. In a lattice formulation one can use the methods of lattice perturbation the- 
ory, which can support the numerical simulations and contain a lot of important information 
concerning the phase structure, masses, renormalized couplings, etc. 

A first look at the phase structure, according to the one-loop lattice effective potential in 
the Landau gauge, shows several qualitative features which are very important from the point 
of view of the existence of the supersymmetric continuum limit: 

• If the quartic couplings Xa,b stabilizing the path integral are not too small, then the 
phase transitions are of second order or very weakly first order. This is the prerequisite 
for the possibility of a continuum limit. 

• The critical value of the ratio r = Xa,b / ' 9 2 above which the first-order phase transition 
spoiling the continuum limit disappears, gets rapidly smaller by decreasing g 2 . This 
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allows us, for small enough values of g 2 , to reach the flow line of renormalization group 
trajectories, which leads to the supersymmetric continuum limit. 

• The wedge containing the AB-phase can be tuned to become narrow, presumably even 
infinitely narrow in the continuum limit. This gives zero mass for the photo-axion, corre- 
sponding to a flat direction in quantum moduli space. The small mass of the photo-axion 
is also supported by the mass estimates from the effective potential (see fig. 7). 

• The broken WT identities corresponding to the global SU(2)k symmetry are very similar 
to the WT identities for flavour non-singlet axial symmetry in QCD with Wilson fermions. 
Therefore, it is plausible that the SU(2)k symmetry is restored in the continuum limit if 
the renormalized fermion mass is tuned to zero. 

These constitute all but one conditions for the existence of the supersymmetric continuum 
limit. The open question is the mass of the photo-dilaton, where the present investigation is 
inconclusive. The left-over combination of the three bare mass parameters, after tuning to a 
zero fermion mass and infinitely narrow AB-wedge, could be a relevant mass parameter for 
the overall length of the vacuum expectation values denoted by v/A = {y\ + v^)? / A. Then 
in the continuum limit the mass of the photo-dilaton would be a function of v/A, which is 
determined by the dynamics and does not need to be identically zero. (In case the width of the 
AB- wedge cannot be tuned to exactly zero in the continuum limit, the same would also hold for 
the mass of the photo-axion.) This expectation seems to be strengthened also by the presence 
of a first-order phase transition at small \a,b-> which can imply a lower limit for the mass of 
the photo-dilaton. Nevertheless, it is very possible that the first-order phase transition and the 
relatively large mass estimate obtained for the photo-dilaton can be avoided by better parameter 
tuning. For instance, one could exploit the above-mentioned infrared attractive direction e 2 in 
eq. (11), which leads to supersymmetric coupling ratios for renormalized couplings and hence 
to a massless photo-dilaton. This question deserves further investigation. 

Further information on the supersymmetric continuum limit can be obtained by numerical 
simulations and by a refinement of the perturbative analysis. After the question of the con- 
tinuum limit has been sufficiently clarified, the predictions of ref. [1] could be exploited in this 
lattice realization of SYM2. The behaviour of the phase transition near the critical point where 
all phases, including the symmetric phase, meet could explain, in the supersymmetric low en- 
ergy effective theory, the absence of a singularity corresponding to the symmetric phase with 
massless elementary fields. The interesting predictions concerning the spectrum of light states, 
including massless monopoles, could be investigated. Since the lattice action contains the pure 
SU(2) Yang- Mills theory as a special case, the understanding of the confinement mechanism in 
the supersymmetric part of the parameter space could be transferred to pure gauge theory. 



Acknowledgements 

I thank Gabriele Veneziano, Zoltan Fodor, Rainer Sommer and Shimon Yankielowicz for helpful 
discussions, and Karl Jansen for correspondence on the effective potential. 



26 



References 

[1] N. Seiberg, E. Witten, Nucl. Phys. B426 (1994) 19; ERRATUM ibid. B430 (1994) 485; 
Nucl. Phys. B431 (1994) 484. 

[2] A. Klemm, W. Lerche, S. Yankielowicz, S. Theisen, Phys. Lett. B344 (1995) 169. 

[3] G. Curci, G. Veneziano, Nucl. Phys. B292 (1987) 555. 

[4] M. Creutz, hep-lat/9411033, to appear in the Proceedings of the Lattice '94 Conference, 
Bielefeld. 

[5] I. Montvay, Phys. Lett. B344 (1995) 176. 

[6] I. Montvay, Nucl. Phys. B (Proc. Suppl.) 26 (1992) 57. 

[7] F. Csikor, Z. Fodor, J. Hein, K. Jansen, A. Jaster, I. Montvay, Phys. Lett. B334 (1994) 
405; 

Z. Fodor, J. Hein, K. Jansen, A. Jaster, I. Montvay, hep-lat/9409017, to appear in Nucl. 
Phys. B. 

[8] G. 't Hooft, Nucl. Phys. B33 (1971) 173; Nucl. Phys. B35 (1971) 167. 

[9] L. Lin, I. Montvay, G. Minister, M. Plagge, H. Wittig, Phys. Lett. B317 (1993) 143. 
[10] M. Bochicchio, L. Maiani, G. Martinelli, G. Rossi, M. Testa, Nucl. Phys. B262 (1985) 331. 
[11] M.F. Sohnius, Phys. Rep. 128 (1985) 39. 



27 



A Appendix: the contribution of the scalar loops 



In this appendix the contribution of the scalar loops to the one-loop effective potential is given 
in the special case where 

= , \ab] = \ab) + 2A,i . (84) 

The notations are somewhat different from those in section 4, namely: 

Xabab = \ab) , AA = (AA) , BB = (BB) , AB = (AB) . (85) 

The expression is ordered according to the powers of 

D A = (m\ + k 2 )- 1 , D B = (m 2 B + A; 2 )" 1 (86) 

and \abab-, ^a- Monomials of these variables are followed in the next one or two lines by their 
coefficients in T> s required in eq. (42): 

1 
1 

Da^abab 
ABB 
D A \ A 
2UAA + 12 BB 

Db^abab 
AAA 

D b \a 
12 AA + 20 BB 

Da 2 ^abab 2 
ABB 2 

Da 2 ^a^abab 
48 AA BB + 32 BB 2 + 16 AB 2 

D A 2 \ A 2 
16 (AA + BB) (7 AA + 3BB) 

DbDa^abab 2 
12 AB 2 -AAABB 

DbDa^a^abab 
80 AA 2 + 64 AB 2 + 32 AA BB + 80 BB 2 
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D B D A \ A 2 
240 (44 + BB) 2 

Db 2 ^abab 2 
AAA 2 

Db 2 ^a^abab 
32 AA 2 + 16 AB 2 + 48 AA BB 

D B 2 \ A 2 
16 (3 44 + 7 BB) (44 + BB) 

D a 3 ^a^abab 2 
16 BB (BB 2 + AABB + 2 AB 2 ) 

Da 3 ^a 2 ^abab 
64 (AA + BB) (2 44 BB + BB 2 + AB 2 ) 

D a 3 \a 3 
64 (3 44 + BB) (44 + BB) 2 

DbDa 2 ^abab 3 
-32 BB (44 BB - AB 2 ) 

DbDa 2 ^a^abab 2 
80 BB 3 + 128 AA AB 2 - 112 AA BB 2 + 224 BB AB 2 

DbDa 2 ^a 2 ^abab 
64 (AA + BB) (7 AA 2 + A AABB + 11 AB 2 + 10 BB 2 ) 

D b D a 2 \a 3 
192 (7 44 + 5 BB) (AA + BB) 2 

Db 2 Da^abab 3 
-32 AA (44 BB - AB 2 ) 

D B 2 D a^A^ABAB 2 

224 AA AB 2 + 128 BB AB 2 - 112 AA 2 BB + 80 AA 3 

D B 2 D a^A 2 ^ABAB 

64 (44 + BB) (lO 44 2 + 4 AA BB + 11 4B 2 + 7 BB 2 ) 

iV^A/ 
192 (5 AA + 7 BB) (44 + BB) 2 

-D^A^A^s^s 2 
16 44 (44 2 + AABB + 2 AB 2 ) 
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Db ^a ^abab 
64 (44 + BB) [AA 2 + 2 44 BB + AB 2 ) 

D B 3 \ A 3 
64 (AA + 3BB)(AA + BB) 2 

DbDa 3 ^abab 4 
— 16 BB 2 [AABB - AB 2 ) 

DbDa 3 ^a^abab 3 
-64 (3 AABB + 2 BB 2 - AB 2 ) (44 BB - AB 2 ) 

DbDa 3 ^a 2 ^abab 2 
64 (AA + BB) [AA 2 BB - 6 AA BB 2 + 5 AA AB 2 + 17 BB AB 2 + 5 BB 3 ) 

DbDa 3 ^a 3 ^abab 
256 (AA + BB) 2 (3 AA 2 + 3 AA BB + 5 BB 2 + 7 AB 2 ) 

D b D a 3 Xa 4 
256 (9 ^4^1 + 5 BB) (AA + BB) 3 

Db 2 Da 2 ^abab 4 
-16 (AA BB - AB 2 ) (3 AA BB - AB 2 ) 

D B 2 D a 2 ^A^ABAB 3 

-64 (44 BB - AB 2 ) (5 AA 2 + 8 AA BB + 5 BB 2 + 2 4B 2 ) 

-Db 2 -D a 2 A 4 2 A 4 s 2 
64 (44 + BB) 2 (7 AA 2 - 20 44 BB + 7 BB 2 + 32 AB 2 ) 

Db 2 Da 2 ^a 3 ^abab 
512 (44 + BB) 2 (7 44 2 + AABB + 7 BB 2 + 11 4B 2 ) 

5376 ( 44 + BB) 4 

Db 3 Da^abab 4 
-16 44 2 (44 BB - AB 2 ) 

D B 3 D a^A^ABAB 3 

-64 (44 BB - AB 2 ) (2 44 2 + 3 44 BB — AB 2 ) 

D B 3 D a^ABAB 2 ^A 2 

64 (44 + BB) (5 44 3 - 6 AA 2 BB + 44 BB 2 + 17 44 4B 2 + 5 BB AB 2 ) 

Db 3 Da^abab ^a 3 
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256 (44 + BB) 2 (5 AA 2 + 3 AA BB + 7 AB 2 + 3 BB 2 ) 

D B 3 D A X A 4 
256 (5 ^4^1 + 9 BB) (AA + BB) 3 
Db 2 Da 3 ^abab 5 
64 BB [A ABB - AB 2 f 

D B 2 D a 3 ^A^ABAB i 

-64 (44 BB - AB 2 ) 

■ (7 AA 2 BB - 3 44 BB 2 - 3 AA AB 2 + 4 BB 3 + 11 BB AB 2 ) 

Db 2 Da 3 ^a 2 ^abab 3 
-256 (3 44 + 5 BB) (AA + BB) 2 (44 BB - AB 2 ) 

Db 2 Da 3 ^a 3 ^abab 2 
256 (AA + BB) 2 

(3 AA 3 - 5 AA 2 BB - 13 AA BB 2 + 18 AA AB 2 + 30 BB AB 2 + 7 BB 3 ) 

Db 2 D a 3 A a 4 ^ 4 BA B 
1024 (44 + BB) 3 (6 44 2 + AABB + 7 BB 2 + 12 4B 2 ) 

Db 2 D a 3 \ 5 
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1024 (9 AA + 7 BB) (44 + BB)" 

64 44 (44 BB - 4B 2 ) 2 

Db 3 D a 2 A 4 A a ba b 4 
-64 (44 BB - 4B 2 ) 

• (4 44 3 - 3 44 2 BB + 7 44 BB 2 + 11 44 4B 2 - 3 BB 4B 2 ) 

Db 3 Da 2 ^a 2 ^abab 3 
-256 (5 44 + 3 BB) (44 + BB) 2 (44 BB - 4B 2 ) 

-Db 3 -D j i 2 A j i 3 A^s^s 2 
256 (44 + BB) 2 

(7 44 3 - 13 44 2 BB - 5 44 BB 2 + 30 4 4 4B 2 + 18 BB 4B 2 + 3 BB 3 ) 

Db 3 Da 2 ^a 4 ^abab 
1024 (44 + BB) 3 (7 44 2 + 44 BB + 6 BB 2 + 12 4B 2 ) 
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1024 (7 AA + 9 BB) (AA + BB) 4 
Db 3 Da 3 ^abab 6 
192 {AABB - AB 2 ) 3 

Db 3 Da 3 ^a^abab 5 
1536 (AA BB - AB 2 ) 3 

Db 3 Da 3 ^a 2 ^abab 4 
-768 (AA + BB) 2 {AA BB - AB 2 ) [AA 2 - AABB + BB 2 + 3 AB 2 ) 

Db 3 Da 3 ^a 3 ^abab 3 
-3072 (AA + BB) 4 {AA BB - AB 2 ) 

Db 3 Da 3 ^a 4 ^abab 2 
3072 (AA + BB) 4 [AA 2 - 3 AA BB + BB 2 + 5 AB 2 ) 

Db 3 Da 3 ^a 5 ^abab 
12288 (AA + BB) 4 [AA 2 + BB 2 + 2 AB 2 ) 

D b 3 D a 3 \a 6 
12288 (AA + BB) 6 
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